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1. Introduction 



Our purpose in this paper is to define and study the properties of a broad generalization 
to n dimensions of the classical Darboux transformation for Sturm-Liouville operators on 
the line. Our approach will stem from a geometric generalization of the basic intertwining 
relations underlying the classical Darboux transformation to the context of certain twisted 
Laplacians acting on the exterior algebra of an oriented Riemannian manifold. 

Let us begin by recalling the essentials of the one-dimensional Darboux transforma- 
tion, [4], [3]. Consider a Sturm-Liouville operator h, given by 

and let be a nowhere vanishing eigenfunction of h with eigenvalue Eq, 

{h - Eq) e-^ = 0. 

The classical Darboux transformation associates to h the Sturm-Liouville operator h de- 
fined by 

^ = -^ + ^(^'- 

where ^ 

V = V + 2x". 

It is straightforward to verify that the operators h and h, shifted by Eq, can be factorized 
in the following way, 

h-EQ = Q+Q-, h-EQ = Q-Q+, 
where and Q~ are first-order differential operators defined by 

ax ax 

The operators H — Eq, H — Eq, and Q~ are therefore related by the intertwining relations 

{h - Eq) g+ = Q+(h- Eq), Q- {h - Eq) = {h - Eq) Q". 

We thus obtain a simple relation between the eigenfunctions (formal or L^) oih and those 
of its Darboux transform h. Indeed, if t/^ is a formal eigenfunction of h with eigenvalue 
E ^ Eq, then it follows immediately from the above intertwining relations that Q~'il; will 
be a formal eigenfunction of h with the same eigenvalue. Conversely, if is a formal 
eigenfunction of h with eigenvalue E ^ Eq, then Q"*"^ will be a formal eigenfunction of h 
with eigenvalue E. It is not difficult to^show that this correspondence also holds at the 
level of the L^ eigenfunctions of h and h, so that^he Darboux transformation establishes 
a correspondence between bound states of h and h. 
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It is well-known that the Darboux transformation plays an important role in the theory 
of soliton solutions of integrable evolution equations and in the method of inverse scatter- 
ing, [5]. It also provides a powerful method for generating new exactly or quasi-exactly 
solvable one-dimensional potentials from known ones [13]. The Darboux transformation 
also appears as a basic tool in the theory of special functions through the factorization 
method of Infeld and Hull, [10]. The problem of extending the Darboux transformation 
to the case of multi-dimensional differential operators is therefore of considerable interest. 

There are at least two natural candidates for what could be called a Darboux trans- 
formation in two dimensions, namely the Laplace, [6], and Moutard, [12], transformations. 
The Moutard transformation is perhaps closer in spirit to the classical Darboux trans- 
formation, since it is based on intertwining relations analogous to the ones given above 
for the Darboux transformation. This will be made explicit in Section 5. On the other 
hand, the Laplace transformation plays a significant role in the realm of integrable sys- 
tems. For example, it naturally gives rise to the Lax representation for the Toda lattice, 
[14]. It is also at the basis of some important recent work on exactly solvable periodic 
two-dimensional Schrodinger operators, [15]. In any case, both the Laplace and Moutard 
transformations preserve the class of linear elliptic second-order differential operators in 
the plane. However, unlike the one-dimensional Darboux transformation, they both suffer 
from the major limitation that they yield only one eigenfunction for the transformed oper- 
ator, [15], the reason being that they do not incorporate the spectral parameter E. Let us 
briefly illustrate this difficulty in the case of the Laplace transformation. The conclusion 
is analogous for the Moutard transformation. Consider the two-dimensional Schrodinger 
equation given by 

in terms of complex coordinates z = + zy), z = ^{x — iy). Under the Laplace trans- 
formation, the wave function i/j gets mapped to the wave function i/j defined by 

7 dip 

and it is easily verified that the transformed wave function '0 satisfies the following Schr5- 
dinger equation for a particle in a magnetic field 



^ d\og{V-E) d ^ 



dzdz dz dz 



= 0. 



The explicit dependence of the coefficient of d/dz on E illustrates our point. There is also 
a drawback which is specific to the Laplace transformation, namely that its geometric gen- 
eralization to n dimensions, [11], applies to a class of highly overdetermined systems which 
bear no relation to any natural spectral problem, although they are of course interesting 
in their own right. Another essential limitation of the Laplace and Moutard transforma- 
tions is that they are only defined for flat Laplacians expressed in Cartesian coordinates, 
whereas many of the Schrodinger operators arising by symmetry reduction involve curved 
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Laplacians in very general coordinate systems. One would therefore like to have a multi- 
dimensional Darboux transformation which is defined in a covariant way, which allows 
for curvature of the underlying manifold and which includes the spectral parameter in a 
natural way. 

An indication on how to proceed is suggested by the work of Andrianov, Borisov and 
loflFe, [1]. In their scheme, one starts from a Schrodinger operator h in the Euclidean plane, 
expressed in Cartesian coordinates, and one constructs, starting from a nowhere vanishing 
eigenfunction of h, the Moutard transform h of h, and a two-by-two matrix Schrodinger 
operator H which splits into the sum of two operators H^^^ and if*^^^. The intertwining 
relations between h, h, H^^^ and H^'^^ imply that h and H^^^ have the same eigenvalues 
except possibly for the zero eigenvalue, and similarly for h and H^'^\ (This construction 
was later generalized in [2] to flat n-dimensional Euclidean space, still using Cartesian 
coordinates in an essential way.) After a careful analysis of this scheme, we conclude that 
the operators h, H and h introduced in [2] can in fact be re-expressed as twisted flat 
Laplacians acting on 0-forms, 1-forms and 2-forms, the latter being identified with 0-forms 
by means of the Hodge operator for the^ underlying two-dimensional flat Euclidean metric. 
The intertwining relations between h, h, H^^^ and H^'^^ will now follow immediately from 
elementary properties of the twisted difl^erentials and codifferentials. 

Starting from this observation, we succeed in our paper to construct a fully covariant 
and coordinate-free multidimensional generalization of the classical Darboux transforma- 
tion, valid on an arbitrary curved n-dimensional oriented Riemannian manifold. Our n- 
dimensional Darboux transformation relates, via intertwining relations involving twisted 
differentials and codifferentials, the spectra and eigenfunctions of a string of n -|- 1 twisted 
Laplacians acting on /c-forms for < k < n. It is noteworthy that these are the twisted 
Laplacians which were used by Witten in his proof of the Morse inequalities based on 
ideas from supersymmetry, [17]. When expressed in any coordinate system, these twisted 
Laplacians take the form of matrix Schrodinger operators acting on the (^) components 
of a A;-form. In the special case where the underlying manifold is flat Euclidean space 
in Cartesian coordinates, our multidimensional Darboux transformation reproduces the 
classical Darboux transformation and the scheme of [1] and [2]. 

Our paper is organized as follows. In Section 2, we define the twisted versions of the 
differential, the codifferential and the Laplacian on forms. When expressed in local coor- 
dinates, the latter will correspond to scalar and matrix Schrodinger operators on curved 
Riemannian manifolds. In Section 3, we first derive the basic intertwining relations defin- 
ing our multidimensional Darboux transformation and give their spectral interpretation. 
This generalization of the Darboux transformation is valid for twisted Laplacians on an ar- 
bitrary n-dimensional oriented Riemannian manifold. We then derive the local coordinate 
expressions of the resulting matrix Schrodinger operators in terms of the seed eigenfunc- 
tion for the original scalar Hamiltonian and the Riemann curvature of the background 
metric. The connection between the spectra of our matrix Hamiltonians admits an inter- 
esting interpretation in terms of supersymmetry, cf. [16], [17], [2], which we briefiy recall. 
In Section 4, we show that our multidimensional Darboux transformation reduces to the 
classical Darboux transformation on the line and to a covariant coordinate-free general- 
ization to curved oriented Riemannian manifolds of the scheme of references [1] and [2] in 
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two dimensions. In Section 5, we derive an n-dimensional generalization of the Moutard 
transformation which apphcs to all twisted Laplacians. The spectral interpretation of the 
multi-dimensional Moutard transformation is significantly more limited than that of the 
multidimensional Darboux transformation, because it only applies to the zero modes of 
the twisted Laplacians. Finally, in Section 6 we obtain new examples of multi-dimensional 
quasi-exactly or exactly solvable matrix Schrodinger operators of physical interest, of which 
very few seem to be known, by applying the multidimensional Darboux transformation to 
various quasi-exactly solvable planar Hamiltonians, [7]. 

2. Twisted Laplacians and Schrodinger operators 

Our purpose in this section is to define a twisted version of the Laplacian on A;-forms on an 
oriented Riemannian manifold. On scalar functions, this twisted Laplacian will correspond 
to a Schrodinger operator. 

We start by setting up some notation. Let M be an n-dimensional oriented Rieman- 
nian manifold, with metric {gij)i<ij<n- We shall denote by f\^{M) the vector space of 
differential /c-forms on M, and by /\{M) — 0^=^ l\^{M) the exterior algebra of differential 
forms on M. In a local coordinate system x = (a;^, . . . , x"^) a A;-form uj e /\^{M) will be 
written as 

^ = l, uj^,...iMd^'" ^■■■^dx''^^ (1) 

where uj^^ is an antisymmetric tensor. Here, and in what follows, we are using the 
summation convention. 

Since M is Riemannian and oriented, it admits a volume form fi e /\^{M), which in 
a positively oriented coordinate system x = {x^, . . . , x^) can be written as 

A* = ^/g{x)dx^ A ... A dx'^, 

where 

^ = det(^^^.). 

The Riemannian metric on M induces a scalar product on the exterior algebra of M in a 
natural way. Namely, if a and /9 are two /c-forms with (antisymmetric) components 0;^^...^^. 
and j^, respectively, one defines 

where the functions 

^ y y ^ji...jk 

are the contravariant components of a, and 
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is the contravariant metric tensor. It is clear that the definition of the function (a, (3) does 
not depend on the local coordinates used; hence one defines the inner product (a, (5) of 
the /c-forms a, /3 e K^{M) by 

{a, (5) = / {a,P) fi. 
Jm 

The inner product {a, (3) will only be defined if the above integral is convergent; for in- 
stance, this will be the case if a and (3 are compactly supported, or in particular if M is a 
compact manifold. 

Given a fc-form 13 e f\^{M), the Hodge star *j3 e l\^~^{M) of f3 is defined as usual as 
the unique (n — A;)-form satisfying 

a ^*|3 = {a,l3) li, \la&f^{M). 

The Hodge star operator * : l\^{M) f\^~'^{M) has the following elementary properties, 
which we list here for future reference: 



(2) 
(3) 



=1, *1 = jJL. (4) 

In local coordinates the components of the Hodge dual of a A;-form a are given by 

(*«W,..,„ = ^S.....«^'^-^^ «eA'(M), (5) 
where e^-^ is antisymmetric under permutations of its indices, and 

ei2...n = +1- 

The Hodge star operator is used to define the codifferential 5 of a /c-form a by the 
formula 

Sa=i-iy^''-^^+Ud{*a), ae/\iM). 

Thus, 6 maps A^(M) into A''"^(M); in particular, 5 = on A°(M) = C~(M). In local 
coordinates, one can show that 

=-V^'«,ii...H_i> (6) 

where V denotes the covariant derivative associated to the Riemannian metric. The odd 
looking sign in the definition of 5 is not arbitrary, since with this choice of sign the operators 
d and 6 are the formal adjoint of one another with respect to the scalar product on A{M), 

4 = 6^+,, 0<k<n-l. 
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It also follows immediately from its definition that 5 is a coboundary operator, namely 

5^ = 0. 

Furthermore, the local exactness of d (Poincare's lemma) implies that 5 is also locally 
exact: in other words, if 5a; = then there is an open neighborhood of every point of M 
on which u = 6a, for some differential form a. 

The Laplacian A : l\{M) A{M), defined by 

-A = d5 + 5d, 

maps each space /\'^{M) into itself for all k = 0,1,..., n. Moreover, from the elemen- 
tary properties of the codifferential it follows that —A is (formally) self-adjoint and non- 
negative: 

(a,-A/3) = (-Aa,/3), a,(3e\\M) 

k 

(cK, -Ack) = {da, da) + {5a, 5a) « ^ A (-^)- 

If / e /\^{M) is a smooth function on M, we obtain using (6) that 

A/ = -M/ = VV./=-L|-(^,..^) (7, 

is the classical Laplace-Beltrami operator on f\^{M). If a; is a /c-form, the components of 
— Aa; in a local coordinate system are given by the following expression: 



k 



q=l 



, to.. ^ . 

»o Jll...lq...lk 



+2 Y. (-ir'-^' -,«....c...c..v w 

l<p<q<fc 

In the above equation a hat over an index means that that index is to be omitted, the 
Riemann tensor R^j^i is defined by 

dvt dri, 

TDi _ ,-ni-nP _ -ni pP 

in terms of the connection coefficients F^^ of the metric, and i?*^ is given in terms of the 
Ricci tensor R^ - by 

TDl TDlh Ah p 

^ j - ^ hj - 9 ^hj- 

A Schrddinger operator (or Hamiltonian) on M is a second-order differential operator 
h : A°(M) ^ A°(M) of the form 

h=-A + V, (9) 



where A is the classical Laplace-Beltrami operator (7), and the potential V is a scalar 
function. Just as the classical Laplace-Beltrami operator is the restriction of —{d6 + Sd) 
to f\^{M), an arbitrary Schrodinger operator h can be expressed in terms of a "twisted" 
version of the latter operator. 

To this end, given a smooth real-valued function x on M we introduce the twisted 
differentials 

and the twisted codifferentials 

The operators d"^ and S"^ have the following properties, which follow directly from analo- 
gous properties of d and S: 

i) (d±)t = 
ii) (d±)2 = (5±)2 = 

iii) d"^ and 5"^ are locally exact 

There are two natural ways of "twisting" the operator A in such a way that the resulting 
operator maps t\'{M) into itself for A; = 0, 1, . . . , n, and is formally self-adjoint and non- 
negative; namely, we can define the twisted Laplacians^ 

H = d-6+ + 6+d- 

and _ 

H = d+6- + 6-d+. 

Let us take the first of these operators, and apply it to a function / e f\^{M). We shall 
need the following lemma, whose proof is straightforward, to simplify the computation: 

Lemma 1. If f and g are scalar functions on M, we have the identity 

5{fdg) = -fAg-Vf-Vg. 

Here V/ denotes the gradient of /, that is the vector field with components 



and V/ • V^f = ((if, dg). Using the lemma we obtain the following expression for the action 
of the operator H on a function / e A°(M): 

Hf = 6+d-f = e^6 [e-^ ■ e-^d{e^f)] 
= e^6{e-^df) + e^6{e-^fdx) 

= -A/ - e^V{e-^) ■ V/ - /Ax - e^Vie'^f) ■ Vx 
= {-A + V-E,)f, 



^ The one-parameter family of twisted Laplacians obtained from H by setting x = 
txo, where t e R and Xo is a Morse function on M , was studied by Witten in his derivation 
of the Morse inequalities using supersymmetric quantum mechanics, cf. [17]. 
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where 

V = (Vx)' - Ax + ^0, (10) 

and Eq is an arbitrary real constant whose purpose will soon become clear. (Note that V 
is real, since x is a real- valued function.) Thus, an arbitrary Schrodinger operator h can 
be represented as 

h^HQ + E^ = {5+d- + d-5+)Q + Eq, 

provided that the function x satisfies (10). The meaning of (10) and of the arbitrary 
constant Eq is apparent if we observe that 

{h - £^o)e~^ = HqC-^ = 6+d-e-^ = 6+ [e-^d{l)] = 0. 

In other words, (10) is equivalent to the fact that is a formal eigenfunction of h 
with eigenvalue Eq. Note that here, and in what follows, by a formal eigenfunction of 
an operator h we simply mean a solution i/j of the eigenvalue equation {h — E)il: — 0, 
regardless of the boundary conditions (like square integrability) that may be used to define 
true eigenfunctions of h. 

Similarly, if we apply the twisted Laplacian H defined above to a scalar function 
/ e A°(M) we obtain 

Hf = 5-d+f={-l^ + V-EQ)f 

with V given by 

V = (Vx)' + Ax + Eq = V + 2Ax, (11) 
since H is obtained from H by replacing x with — x- Letting 

h=-A + V (12) 

we then have 

{h - EQ)e^ = 0. 

Note that the Schrodinger operator h is the Moutard transform of /i, [12] . 
3. The Multidimensional Darboux Transformation 

The goal of this section is to define the multidimensional Darboux transformation. The 
definition will be naturally suggested by some fundamental intertwining relations involving 
the differential operators on /\{M) introduced in the previous section. 

Let us consider the operator H : /\{M) — > /\{M), which we shall decompose as follows: 
where 

H^'^=d-5+, H^^^ = 5+d-. 
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Observe that, since S'^ and d are coboundary operators, we have 

^(1) ^(2) ^ ^(2) ^(1) ^ Q_ (13) 

Since the operator H maps f\^{M) into itself for A; = 0, 1, . . . , n, we can also write 

n 
fe=0 

where is the restriction of H to f\^{M). The operator can be decomposed as 

where and -ff^^'* are the restrictions of H^^^ and iJ*^^^ to /\^{M)^ namely 

o-(i) _ J- x+ o-(2) _ r+ j- 

Notice that 



= = 0. 



Prom (13) we obtain 



= = 0. (14) 



By construction, H is formally self-adjoint and non-negative. Therefore the same is true 
for the operators : A''{M) A''{M) for A; = 0, 1, . . . , n. Likewise, both hI^^ and i^f ^ 
are formally self-adjoint and non-negative, the latter property being a consequence of the 
identities 

{a,Hj^^a) = {d~a,d~a), 

{a,Hl^^a) = {S+a,6+a), ae/\{M), k = 0,l,...,n. 

The following intertwining relations are an immediate consequence of the definition 
of and if (2). 

ffW d- = d-H^'\ 

^Ui <li = H^k^ ^Ui (15) 
HiZ dk = d^H^k\ = 0, 1, . . . , n - 1. (16) 

The intertwining relations (15)-(16) have important consequences for the spectra of the 
operators h'^^.^ and h'^^ , that we shall now explore. 
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These imply 



Proposition 2. 

i) If u) e l\^'^^{M) is an eigenform of H^'^.^ with, eigenvalue X 0, then = (^^+i<^ is 

(2) 

an eigenform of with the same eigenvalue. 

ii) Likewise, ifuJE/\ (M) is an eigenform ofHf, with eigenvalue X ^ 0, then d~u> = df^uj 
is an eigenform of H^^i with the same eigenvalue. 

Proof: The only point that is not an immediate consequence of the intertwining relations 
is that in part i) cannot vanish identically, and likewise for d~uj in part ii). Let us 
show, for instance, that 5^uj 7^ in part i). If ^"'"0; = 0, then by the local exactness of 5^ 
(property iii) in Section 2) for every p G M there is an open neighborhood of p such 
that uj = S'^ap on Up, for some differential form on Up. But this would imply that 

-^fe+i'^ = d~S~^{S~^ap) = on Up, for every p e M, so that i^^^^o; = on M. Hence 
A = 0, contradicting the hypothesis. Q.E.D. 



Let L be a differential operator on l\{M) mapping each subspace f\^{M) into itself. 
Examples of such an operator are H, H^^^ and H^'^K We shall denote by A{L) C l\{M) 
the set of admissible forms for L, that is the subspace of the space of square- integr able 
forms L^(^/\{M)^ satisfying any additional boundary or asymptotic conditions that are 
appropriate for the problem being considered. The set of admissible forms for is then 
^(Lfc) = ^fc(L) = A{L) n K^{M). The spectrum of L, denoted by cr(L), is the set of 
numbers A such that there is an admissible non-zero eigenform ui G A{L) satisfying the 
eigenvalue equation Luj = Xuj. We shall also use the convenient notation 

a'{L) = a{L) - {0}. 

The linear space of all /c-forms uj satisfying the eigenvalue equation Luj = Xuj will be 
denoted by Aa(-^)- Equivalently, Aa(-^) is the set of eigenforms of with eigenvalue A, 
together with the zero fc-form. 

By the formal self-adjointness and non-negativity of Hj^, h'^^ and h'^\ the eigenval- 
ues of these operators are real and non-negative. Proposition 2 has the following immediate 
corollary: 

Corollary 3. If the operators 5^ and d~ map A = A{H) into itself, then the spectra of 
nj^^ and h'^^i are related by 

a'{HP) = a'iHi'^,), k = 0,1, . . . ,n - 1. (17) 

In the remainder of this section, we shall not distinguish between true and formal 
eigenforms, unless otherwise indicated. The following lemma, whose proof is elementary, 
will have non-trivial consequences in what follows: 
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Lemma 4. Let V be a vector space, and let L : V ^ V be a linear operator. Suppose 
that L = + L2, with L2 = L2 = 0. The following statements are then true: 

i) If V is an eigenvector of L with eigenvalue X, either L-v = Xv and L-v = for some 
i, j G {1,2} with i j, or L^v is an eigenvector of with eigenvalue X for i = 1,2. 

ii) a{L) Ca{L^)Ua{L2). 
Hi) a'{L) = a'{L,)Ua'{L2). 

From the previous lemma and (14), it follows that 

a'{H,) = a'{Hi'^)Ua'{HP). (18) 

Thus, the spectrum of iJ^, with the possible exception of the zero eigenvalue, is simply 

the union of the spectra of its components Hj^^ and The spectra of Hj^^-^^ and Hj^^ 

(k = 0, 1, . . . , n — 1) are identical, except perhaps for the zero eigenvalue. The operator 

S'^ maps eigenforms of ffj^^-^ with non-zero eigenvalue into eigenforms of Hj^^ with the 

same eigenvalue, and similarly d~ maps eigenforms of with non-zero eigenvalue into 

eigenforms of Hj^_^.^ with the same eigenvalue. From the identity (17) we easily obtain the 
following relation between the even and odd components of H: 




Definition 5. Let A 7^ 0, and /c = 0, 1, . . . , n - 1. If a; e aJ^^(^^'^^^), its Darboux 
transform is the A;- form S'^uj E /\^(^H^'^^y Similarly, the Darboux transform of a A;-form 
OJ e Aa(^^^^) is the (A; + l)-form d-co e A'1+^{H^^^). 

From the definition of H^'^^ it directly follows that if A ^ the restriction 5+ : 
yyfc+i (ij (1)) /\^ (iJ (2)) is invertible, its inverse being the restriction A~^ d~ : Aa {H^^'>) 

The Darboux transformation we have just defined acts in a natural way on eigenforms 
of the partial Hamiltonians H^^^ and H^'^\ In the same spirit, we shall see next how to 
use the Darboux transformation to construct new eigenforms of H of degree k — 1 and/or 
k + 1 starting from a given eigenform of degree k with non-zero eigenvalue. To this end, let 
u) G l\^{M) be an eigenform of H with eigenvalue A 7^ 0. By Lemma 4, either H^'^^ u> — Xlu 
for some i G {1, 2}, or uj is an eigenform of i?*^*^ with eigenvalue A 7^ for z = 1, 2. In 
the first case, the Darboux transform of to is well defined according to Definition 5, and 
belongs to either AA~^(-f^^^^) (when i = 1), or to Aa"^^ (-ff^"^-*) (when i = 2). By Lemma 
4, it follows that the Darboux transform of u is an eigenform of H in this case. In the 
second case, the Darboux transforms of both H^^^u) and H^'^^u) are defined and, as before, 
are eigenforms of H with eigenvalue A. Hence in this case both oj = 5'^d~5'^ oj and 

d~H^'^^oj = d~5'^d~ OJ are eigenforms of H with eigenvalue A and degree equal to A; — 1 
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and + 1, respectively. Note that in this case both H^^^u and H^^^co are also eigenforms 
of H of degree k with eigenvalue A, neither of which is proportional to u (although the 
span of H^^'^u and if^^) 

u is obviously two-dimensional). Thus in the second case, i.e, 
when H^^^ a; ^ for i = 1, 2, we can construct three new eigenforms of H of degrees k — 1, 
k and k+1 and eigenvalue A starting from a known eigenform a; G Aa {H) . 

The above construction could have been carried out using the twisted Laplacian H 
instead of H. However, we shall now show that the two constructions are equivalent: 

Proposition 6. For k = 0,1, . . . ,n, the operators H^_^ and H^^ are linearly equivalent 
under Hodge duality: ^ 

Hk = i*r'H^-k*- (19) 
Proof: An elementary calculation shows that 

from which (19) follows using (3). Q.E.D. 

In other words, if we identify the space l\^~^{M) to l\^{M) under the Hodge star oper- 
ator, H^_j, is transformed into Hj^. In particular, is equivalent under this identification 
to Hq = h — Eq, the Moutard transform oi Hq = h — Eq. 

We end this section by deriving local coordinate expressions for thejnultidimensional 
Darboux transformation and for the component Hamiltonians and i?^. By definition, 

d~uj — e~^d(e^uj) — dw -\- dx f\ oj , 

so that if a; e A^{M) we have the following local coordinate expression 

k 
r=l 

Here the square brackets denote antisymmetrization, and Xj = ^jX = dx/dx^ . Similarly, 
using (6) we obtain 

(5+'^)n....... = (-V^' + x^)^,n.....-.' A'W- (21) 

Using the local formulas (20) and (21) we obtain, after a straightforward calculation 
iH,u;),^„, = [-V^.V^- + (Vx)^ - Ax] u;,^..,, 

k 

+ 5 ^ir) X ' ^ii...ir—ljir+l---ik 

r=l 

13 



k 

This formula can be further simphfied as foUows. In the first place, note that since x is a 
scalar function we have 

{ V, J X = 2 V,^x = 2 V,^ V^x- 
Secondly, from the identity 

k 
r=l 

we obtain, after straightforward manipulations 

k k 

l<r<q<k 

Substituting the previous formulas into (22) we get the following coordinate expression for 
the action of i?^: 

(^.'^)n...u= + (VX)^ - Ax] c.,,...,, 

fc 

.t;.r..,.. 

r=l 

+ E (-i)^^^+^^^.\-..,...c...?;....,- (23) 

l<r<g<fe 

When X = the previous expression reduces to the formula for minus the Laplacian, 
eq. (8). In particular, the operator Hj^. has the structure 

Hk = -A, + V„ 

where is the restriction of the Laplacian to l\^{M), and Vjr. acts like a matrix potential 
on the components of any A;-form: 

k 

{yk^)i,...iu = [(Vx)' - Ax] cj,^...,^ +2^ ^'^i^X-^i,...ir-,3i.+.:.iu 

r=l 
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k 
r=l 

The local coordinate expressions for the component Hamiltonians Hj. are obtained from 
the previous formulas by replacing x with — 

As first pointed out by Witten, [16] , [17] , [2] , the connection between the spectra of the 
component Hamiltonians Hj^ discussed in this section admits an interesting interpretation 
in terms of supersymmetry. Indeed, the n + 1 homogeneous components a;^ of degree k of 
a differential form Vt = 0^=0 '^k interpreted as the components of a supermultiplet, 

with /c-forms regarded as being bosonic or fermionic depending on whether k is even or 
odd, respectively. The supercharges are by definition the operators 

Q- = 5+, Q+ = {Q-)^ = d-, 

while the supersymmetric Hamiltonian 

H^{Q+,Q-} 

is just the twisted Laplacian. The remaining commutation relations defining the standard 
supersymmetry algebra 

{g±,Q±} = [g±,iy] = 

hold thanks to the elementary properties of d^, S^. 

We can also introduce fermion creation and annihilation operators 

where _l denotes the inner product. The usual fermionic anticommutation relations 

{b^+,lri+} = {b-,b-} = 0, {b'+,b-} = S], 

as well as the identity 

b'+ = (bn^ ^i9''b-y. 

follow easily from well known exterior algebra identities. The supercharges can be 
expressed in terms of the creation and annihilation operators as 

where 

Qf^ = y_ (tV^ + V^x) • dx'' A ... A dx''^ = fify. q^^uj. 

One can easily check the identity 

H=-A + {Vxf + V,V,.x • [b'+,bi-], 

which generalizes formula (13) of [17] (where an orthonormal basis of the tangent space is 
used to define creation and annihilation operators). 
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4. The Darboux Transformation in Low Dimensions 



Let us begin by showing that in the one-dimensional case the multidimensional Darboiix 
transformation reduces to the classical Darbouxjransformation. Indeed, in this case M = 
IR, Hq = h — Eq, and is equivalent to Hq = h — Eq under the Hodge duality. If 



then from (11) we have 

where 
and 



V = V + 2x" 



{h- EQ)e-^ = 0. (24) 



The operator d maps eigenfunctions of h = Hq + Eq with eigenvalue E ^ Eq into 
eigenfunctions of H-^ with eigenvalue E — Eq. Therefore, the operator Q~ — *d~ will 
map eigenfunctions of h with eigenvalue E ^ Eq into eigenfunctions of h with the same 
eigenvalue. If a; is a cartesian coordinate we easily obtain 

SO that 

Q- = ^+x'. (25) 

Similarly, if (j) is an eigenfunction of h = H + Eq with eigenvalue E then *0 will be an 
eigenfunction of with eigenvalue E — Eq, so that 5+(*0) is an eigenfunction of h with 
eigenvalue E. Proceeding as before we have 

<^+(*<^) = (-^+x')0, 



so that the operator 



Q- = ~W (26) 



maps eigenfunctions of h with eigenvalue E ^ Eq into eigenfunctions of h with the same 
eigenvalue. From the properties of d"^ and 5^, it follows that Q~ and are formally the 
adjoint of one another under the Euclidean scalar product on R, and that 

h-EQ^Q+Q-, h-EQ^Q-Q+. (27) 

Equations (24), (25), (26) and (27) express the classical Darboux transformation. 
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The two-dimensional Darboux transformation generalizes to curved oriented surfaces 
the Darboux transformation introduced in [1] for in cartesian coordinates. Indeed, let 
M be a two-dimensional oriented Riemannian manifold. The component Hamiltonians 
are in this case Hq ^ h - Eq, : /\^{M)^^ A^(M) and : A^{M) A^(M), which 
is equivalent under the Hodge duality to Hq = h — Eq. The scalar Hamiltonians h and 
h are given by eqs. (9), (10), (12) and (11), while the function as usual satisfies 
{h — Eq) = 0. To the operator acting on one-forms there corresponds an operator 
Hi acting on vector fields, defined by 

H^X = (iyiX^)«, (28) 

where = g^jX^dx'^ is the one-form associated to the vector field X = X'^d/dx'^, and 
tl = \?~^. Using (23) we easily find the following local coordinate expression for H^X: 

{Hixy = [-V,. V^' + (Vx)' - Ax] X' + 2 V^V^-x • x^ + i?} x^. 
Prom the well-known two-dimensional identity 

where K is the Gaussian curvature of M, we obtain the equivalent expression 

{HiXy ^ [-V^.V^' + (Vx)' -Ax + K] X' + 2 V^V,.x • X^. (29) 

In flat space and cartesian coordinates, the above formula for reduces to formula (9) of 
[1]. By Corollary 3 and eq. (18), the spectra of the Hamiltonians h, h and Hi are related 

by 

a\Hi) = a'{h - E^) U a'{h - E^). 

We shall now derive an expression for the two-dimensional Darboux transformation 
in local coordinates. Let the operators qf and be defined in local coordinates by 



<lf = TV, + X, = Te^^V.e^^ 

pf = Vge^jq'^- (30) 

In particular, notice that the operators (resp. pf) transform like the components of a 
covariant tensor (resp. pseudo-tensor) of rank one under changes of local coordinates. 

If '0 is an eigenfunction of h with eigenvalue E ^ Eq, then d~'il^ is an eigenfunction of 
Hi with eigenvalue E — Eq. Using the general formula (20) we have 

(d-^), = q-^. (31) 
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From (28) it follows that the vector field with components 

g'^qT^I; = g^-V; (32) 

is an eigenvector of with eigenvalue E — Eg. Suppose now that (f) is an eigenfunction of 
h with eigenvalue E ^ Eq. Then *0 is an eigenform of H2 with eigenvalue E — Eq, so that 
5+ (*</)) is an eigenform of with eigenvalue E — Eq. From the local formulas (5) and 
(21) and the fact that the pseudo-tensor with components ^/ge^j is covariantly constant, 
it follows that 

5+ {*(/)) i = 

Therefore the vector field with components 



g^^pj<j> = p'-<j> (33) 
is an eigenvector of with eigenvalue E — Eq. Note also that can be expressed as 



Conversely, if = dx^ is an eigenform of H[ ' with eigenvalue A 7^ 0, then (Lemma 
4) \1/ is an eigenform of with the same eigenvalue, and is an eigenfunction of h 
with eigenvalue X + Eq. From (6) we obtain 

= -e^v'ie-^iP,) = q'+iP, = q+iP\ (34) 

where the vector field = ip^d/dx^ is an eigenvector of H^. Likewise, if ^ is an eigenform 

(2) 

of Hi with eigenvalue A 7^ then \1/ is an eigenform of with the same eigenvalue, 
and (J+^f is an eigenfunction of H2 with eigenvalue A + Eq. It follows that is an 

eigenfunction of h with eigenvalue X + Eq. Using the local coordinate formulas (5) and 
(20) we easily obtain 

*d-^r = p+^\ (35) 

where again the functions ip'^ are the components of an eigenvector of Hi with eigenvalue 
A. As before, in flat space and cartesian coordinates eqs. (30), (32)-(35) reduce to the 
corresponding formulas in [1] . 

The above formulas expressing the two-dimensional Darboirx transformation in terms 
of the operators and pf suggest that the component Hamiltonians and the intertwining 
relations can also be written in terms of the latter operators. A straightforward compu- 
tation using the local coordinate expressions for *, d and 6 shows that this is indeed the 
case. More precisely, we have 

h-EQ = q+q'- , h-EQ = q-q'+ = p+p'' 
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and = h[^^ + h[^\ with 



3 

\ i 



Notice that, strictly speaking, [H^^'>Y ■ is not the (z, j)-th matrix element of the operator 

, since V ^il^^ in general depends on all the components of ^. Similarly, the intertwining 
relations (15)-(16) can be written as 

4{H^'^)\ = {h-E,)q-^ 
{HWy.q^- =q^-{h-E,) 
pnH^'^yj = ih-E,)p+ 

In flat space and cartesian coordinates, the above expressions reduce to the corresponding 
ones in [1]. 

We shall now show how the classical Moutard transform, [12], is generalized in the case 
of an oriented Riemannian surface. To this end, suppose that t/j is a, formal eigenfunction 
of h with the same eigenvalue Eq as e~^. Note that ^/^ need not be proportional to e~^, 
since we are dealing with formal eigenfunctions. A Moutard transform of •0 is any function 
ip satisfying 

d+'4j = S-{*i(;). (36) 

Locally, the above equation has a solution if and only if [h — Eq)iI) = 0. Indeed, by the 
local exactness of d'^ the compatibility condition of (36) is 

d+5-{*ij) = 0. 

The Hodge dual of the latter equation is simply 

S+d'ip ={h- Eq)iP = 0, 

as claimed. Note that the Moutard transform ip is locally defined by (36) only up to a 
constant multiple of e^. Indeed, if t/^j^ and 1(^2 Moutard transforms of then their 

difference satisfies ^"''('01 — ■^^2) = 0) that is d e~'^{^-^ — ■!^2)j = 0. 
If '0 is any Moutard transform of 1/', then 

(h - E^)^ = 0, 

where as before h is the Moutard transform of h. Indeed, applying 6~ to (36) we obtain 

= 6-d+{^) = Hq4; ^{h- Eq)^;. 
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5. The Multidimensional Moutard Transformation 

We shall derive in this section a generalization to oriented Riemannian manifolds of ar- 
bitrary dimension of the classical two-dimensional Moutard transformation introduced in 
the previous section. The key to this generalization is a remarkable connection between 
the zero eigenspaces of the operators and H^^^ (or, equivalently, of the components 
H^^^ and H^j^l^ of H defined below) that we shall describe next. 

Throughout this section, M will denote an oriented Riemannian manifold of dimension 
n > 2 with trivial de Rham cohomology. In particular, the latter condition will always 
hold if we restrict ourselves to a contractible open subset of M. Let us decompose H as 

with ^ ^ 

so that _ 

i/W = z = l,2, 

as in Proposition 6. As before, the operators H^'^^ are formally non-negative and are the 
formal adjoint of one another. If A; = 2, 3, . . . , n, let w e Aq (^^^^) be a zero mode of ' , 
i.e, a fc-form u satisfying the equation H^-^'^u = 0. In this section we shall exclusively deal 
with formal eigenforms, so that in particular u> is not required to be square-integrable. By 
analogy with (36), it is natural to consider the equation 

d+u; = 5-u; (37) 

as an equation for the {k — 2)-form U. The properties of (37), which are easily established, 
are analogous to those of (36), namely: 

i) Equation (37) is compatible if and only if w is a zero mode of H^'^\ 

ii) lJ is uniquely defined by (37) up to an element of lm.d^_^ — Im(i^_3. In particular, 
if A; = 2 we have Im d_i = R (the space of constant functions on M), so that Im d'^^^ = 

iii) Any solution U of (37) is a zero mode of -^^^^2' have 

From the above properties it follows that (37) defines a mapping 

which is easily seen to be an isomorphism by the assumption on the de Rham cohomology of 
M. Note that for the above mapping to be non-trivial M must be non-compact. Indeed, if 
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M is compact then an easy integration by parts argument implies that Ao(-f^^^^) = Ini<^, 



_ k+l 

and A5"'(i?^')) =Ini4_3. 

Totally analogous considerations can be made for the components Hj^^ and Hj^_^2 of 
H. To be precise, consider the equation 

5+u = d-u;, (38) 

where a; is a A;- form and A; = 0, 1, . . . , n — 2. As before, the integrability condition for (38) 

is that a; be a zero mode of Hj^\ in which case ui e /\^'^'^{M) is uniquely defined modulo 

Im5^^3 = e^Im5;^^.3 (with \m.5^_^i = *{Im6_i) = M//), and is a zero mode of Hj^^^- 
Equation (38) thus defines a mapping 

which is again an isomorphism. Clearly, the mappings ~ and " are related by Hodge 
duality. More precisely, an elementary calculation yields the following result: 

Proposition 7. The maps (")^ and are conjugated under Hodge duality: 

C)k = -i*r'<-)n-k°*- (39) 

Comparing equations (36) and (37), we see that a Moutard transform of a function 
ijj G Ao(-f^''^'*) = Ao{H) on a two-dimensional manifold M is simply any function ijj in the 

equivalence class of which will be an element of Ao(-f^*^^'*) = /\o{H). This observation 
motivates the following general definition: 

Definition 8. The Moutard transform of w G /\q {^H ^'^^^ / Im d'i^_^ is the element 

-k-2. 



ofKr'-\H^^))/i^du_,. 

The domain and the range of the Moutard operator defined above follow immediately 
from the identity 

or equivalently, using (39), 

~ = -(*o-). 

In other words, if a; G Aq I'^'^fc-i its Moutard transform is the unique solution 

~ g A''~^~^(M)/Im(i+_^_3 of the equation 

d+uj^5-{*uj), (40) 

which is automatically an element of S^~^~^ {H^^^^ / lm.d'^_y,_^. As before, for the Mou- 
tard operator to be non-trivial M must be non-compact, since otherwise Ao(-f^*^^^) = 

Im(i^_j^ and Ao~'^~^ (-^'^^'*) = ^^d'^_^_^. Notice also that the Moutard transform of a k- 
form uj has different degree than cu, unless n = 2m is even and k = m — 1, with m = 1, 2, . . .. 
For instance, when m = 1 we have n — 2, k — n — k — 2 — 0, and we obtain the general- 
ization of the classical Moutard transformation to Riemannian surfaces introduced in the 
previous section. 
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6. Examples 

We present in this section a few examples of the two-dimensional Darboiix transformation 
on curved surfaces based on the theory of quasi-exact ly solvable Hamiltonians, [7], [8]. 

As our first example, consider the first-order differential operators 

The above differential operators span a Lie algebra g ^ 0t(2, M) k M^, which coincides with 
the canonical form 1.11 in the classification of Lie algebras of differential operators in two 
variables of reference [9] for n = r = 2. The latter Lie algebra is quasi-exactly solvable, 
[7], since it preserves the subspace jV C C°°(M^) whose elements are the polynomials in 
the variables {x, y) of total degree less than or equal to 2. 

Let J denote the following element of the universal enveloping algebra of g: 

J = ( Ji)2 + ( j2)2 + 4 (a - 6) ( J3)2 + 4 5 ( J5)2 + 4 « _ 26) {J^f 

+ 2 a { J^ J^} - 4 (2a + 6) + 2 (a - 126) + 5(76 - a) , 

where a, 6 are real parameters such that 

a > 26 > 0. (41) 

If we define 

(7= ^log[l + 2(a-26)a;2] - | log[l + 2(aa;2 + 26^^)] , 
then it can be shown that 

-e" -J ■ e-" = h, (42) 

where h — —A + V{x,y) is a Schrodinger operator on the manifold M — M.'^ endowed with 
an appropriate metric. More precisely (cf. [7], Example 4.3.3), the contravariant metric 
tensor {g'^^) has components 

= (1 + 2 a x^) [1 + 2{a - 2b)x^] , 
g^"^ = 2axy[l + 2{a- 2b)x^], 
g^^ = l + 46y2 + 4a(a-26)a;2y^ 

and Gaussian curvature 

K=-2a [l + 4(a-26) x^)], 
and the potential V{x,y) is given by 

6 48 6 



V = -3a{a-2b)x' 



l + 2(a-26)a;2 1 + 2 (aa;^ + 26y2) ' 
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Since J belongs to the universal enveloping algebra of Q by construction, it restricts to the 
finite-dimensional vector space jV. We can therefore easily diagonalize JjA/", obtaining the 
eigenvalues 

Aq = -5a + 236 + 4s, = -5 a + 23 6 - 4s, A2 = -a + 36, 

A3 = -9a + 276, A4 = -3a + 156, X^ = -3a + 7b, 

where 

s = \/o2 -2a6 + 962. 
Their corresponding eigenfunctions are given by 

<^Q — 3b + s + 2a {a — 2b) x^, ^p-^ = 3b — s + 2 a {a — 2b) x^, 
(P2 = l + 2ax'^ -16by'^, (p^ = x, <P4 = y, (p^ = xy. 

By (42), the Schrodinger operator h possesses the six eigenfunctions 

iPi = e''^,, 0<z<5, 

with energies 

E, = -\. (43) 

The first of these eigenfunctions has no zeros, and therefore it must correspond to the 
ground state of h. The ground state energy of h is then given by 

EQ = ba-23b-4s, 

which is indeed manifestly lower than the remaining five energies (43) by (41). 

To define the twisted Hamiltonian H and the Darboux transform of the eigenfunctions 
V'j, it is convenient to take as a constant multiple of the ground state i/jq, since this 
eigenfunction has no zeros. We shall therefore define 

X = - log(a;^ + c) - cr, 

where the constant c > is given by 

3b + s 
^~ 2 a {a -2b) ' 

According to (29), the action oi H on a. vector field X can be expressed as 

{H^xy = i-VjV^ + V^) X' + V'j X\ 
23 



where = V — Eq + K and V^j = 2 V* V^- x- After a long but straightforward calculation 
we find that 

= -7a + 23b + 4:3-11 a-fx'^ - ^ ^ - ^, ^, „, , 



48abxy{l + 2-fx'^) 
' l + 2{ax^ + 2by^) ' 
2 4:8axy {b + tt'j x'^) 



a;2/Pg(a;) 



where 



and 



1 + 2(ax2 + 26y2) 7 a (c + a;2)2 (1 + 2 7^2) 
7 = a - 26 > 



Pg (x) = 48 7^ a;^ + 4 a 7^ (9 a + 12 6 + 4 s) + 4 7 (5 - 9 a 6 - a s + 18 6 s + 54 r ) 
+ 5a^-22a6-4as + 30 6s + 90 6^ 



Note, in particular, that ^^i- 

The Darboux transforms d~'i(j- of the eigenfunctions t/j-, 1 < i < 5, constructed above 
can also be computed in a straightforward way using (31), with the following result: 



4 e'^s X 



c + a;2 



dx 



d ■02 = 2 e 
d~^3 



x{16by'^ + 2ac-l) 



c + a;2 



dx — 16 bydy 



d V4 
d-ip^ 



^'—^dx 
c + x^ 

2xy 



= e 



c + a;2 

y{c-x'^) 
c + a;2 



dx + dy 
dx + xdy 



As explained in Section 4, the five vector fields associated to these one-forms are formal 
eigenvectors of with eigenvalues E- — Eq, 1 < i < 5. For example, the first of these 
vector fields is ^ 

4sa; ^^^\^ e" [(1 + 2 ax^) 9^ + 2 a x 2/9 1 . 

In this case, it is easy to show that the five formal eigenvectors of Hi constructed above 
are actually square-integrable. Indeed, the square norm of a vector field X can be 
expressed as 

1^11'= / ^g'^X.X^dxdy, 
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where g = det(^^j) and [i = 1,2) are the components of the associated one-form . 
But 

e^" ^= (1 ^2ax^ + 4by'^)~'^^'^ < const. r~'^, r = x'^ + y'^ ^ oo, 

and a straightforward calculation shows that e~'^'^ g'^^ X^X^ is bounded by a constant times 
r'^ as r ^ oo when X'^ = d~'ip^ for i = 1, . . . , 5, thus proving our contention. 
The limiting case 

a = 2b 

of the previous example is worth studying, because in this case the curvature is constant 
and negative: 

K = -46 < 0. 

Diagonalizing again J\J\f we find the following eigenvalues 

Ao = 25 6, Ai = 9 6, A2 = 6, 

of respective multiplicities 1, 2 and 3. The corresponding eigenfunctions are 

<^o = l' ^i=x, (P2 = y, ^3 = xy, (p^ = 12bx^ -1, (p^ = 12by^ -1. 

As before, the six functions i/j^ = e'^ (f^ are eigenfunctions of h with eigenvalue = — A^, 
where cr can be taken as 

3 

a = --log(l + 46r^). 

Again, the first of these eigenfunctions has no zeros, and therefore corresponds to the 
ground state, with ground energy 

= -25 b. 

Taking i/jq = e~^, or equivalently 

X= -cr- 

we find 

^1-^^^ l + 46r2' ^ 2-V 1- 1 + 4^,^2- 

The Darboux transforms of the eigenfunctions i/j^ (1 < * < 5) are easily computed. Their 
associated vector fields, which are as usual eigenvectors of with eigenvalues E^ — Eq, 
are 

e'^ [{l+4:bx^)d^ + 4:bxydy\, [Uxyd^ + {1 + Uy^) d^] , 

[(l + 86a;2)ya^ + (l + 86/)a;aJ , 24be'' {{\^ 4bx'')xd^^ 4bx^ yd^ , 
246e^ [46xy2a^ + (l + 46y2)^5j . 

Since 
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an argument analogous to the one for the previous example shows that the latter eigen- 
vectors are all square integrable. 

Consider now the canonical form 2.3 (with n = 2) in the classification of reference [9], 
which is the Lie algebra Q ~ sl{3, M) spanned by the first-order differential operators 

J'^ = x^d^ + xydy-2x, = xy + dy - 2y . 

The latter Lie algebra leaves invariant the finite- dimensional subspace J\f C C°°(M^) in- 
troduced in the previous examples, so that it is again quasi-exactly solvable. Let J denote 
the element of the universal enveloping algebra of g given by 

J = a [{J'f + e {J'f] + c [e {jy + {J^] + b {{J\ J'} + { J2, J8} + 7 J3 + 7 J') , 

where a, c, £ are positive real parameters, 6 > 0, and 

7 = £ac-6^>0. (44) 

If 

p = b + c{ex'^ + y^) (45) 

and 

5 7 
= - g log(P^ +7) + 4 ^7~^^^ arctan(7"^/2p) , 

it is shown in [7] that (42) defines again a Schrodinger operator h, with potential V given 

by 

71 , 3 (192b'^ -45aci) p+U3abc£-192b^ 
V = — — b — - p + 



4 4 4 (p2 + 7) 

The metric of the Riemannian manifold M = has now contravariant components 

g^^ = a + x^ (b + p), 

g^^ = ea + y''ib + p), 

and Gaussian curvature given by 

K = -2p. 

The eigenvalues of J\J\f are easily found to be 

Ao = 46 + 2s, Ai=46-2s, A2 = 126, A3 = 26, 

where 

s = Vl6b'^ + 2ac£, 
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and the last two eigenvalues have multiplicity 2. The corresponding eigenfunctions are 
respectively 

c c 

On account of (42), multiplying these eigenfunctions of J by the factor e"' we obtain six 
eigenfunctions '0^ = ip^ of h, with energies E- = — A^. As in the previous examples, 
the first of these eigenfunctions never vanishes, and therefore it corresponds to the ground 
state of h, with ground energy given by 

Eq = -4b -2s. 

We therefore take 

X = -logV'o = -(J -log{£x'^ + y'^ + k), 

with 



obtaining 



k = > 0, 



„ 55, 11 (1926^ -45ac£)p+ 143a6c£- 1926^ 
V^ = 2s-—b-—p + 



4 4 4 (p2 + ^) 

and 

1 2cxyP^{p) 2 _fyi 

' (p2 + 7)(p + s-56)2' 1 ''^2' 

where 

Pg (p) = (27 6 - 4 s) p3 + (9 a c£ + 14 6 s - 79 6^) 

+ (3176^-766^5- 16a 6c£+ 6a c£s)p 
-2016^ + 187a62c£- 5067s + 14a^ . 
Note that the denominator of V^2 never vanishes on account of (44)-(45). 

The Darboux transforms of the eigenfunctions ip^ {I < i < 5) are given by 

4 s e*^ 

= /, , . 2 I 2^ i^xdx + ydy), 
c[k + £x'^ + y^) 

^"^2 = , , .^""s I 2 [^2;(/c + 2|/2)(/a;-(/c + 2^a;2)2/d2/] , 

^"V^4= ;,^^^2^^2 [{k-£x^ + y^)dx-2xydy] , 
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The vector fields associated to these one-forms are eigenvectors of with eigenvalue 
E^ — Eq. As in the previous examples, the latter vector fields are easily found to be square 
integrable on M on account of the asymptotic behavior at infinity of e^'^ ^ and e~'^ [dT'^^. 
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